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Abstract 

The  idea  of  replacing  a  divergence  constraint  by  a  divergence  boundary 
condition  is  investigated.  The  connections  between  the  fonnulationfl  are 
considered  in  detail.  It  is  shown  that  the  most  common  methods  of  using 
divergence  boundary  conditions  do  not  always  work  properly.  Necessary 
and  sufficient  conditions  for  the  equivalence  of  the  formulations  are  given. 
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1  Introduction 

A  well  known  problem  in  computational  electromagnetics  is  the  appearance  of 
‘Spurious  modes”  arising  from  incorrectly  imposing  divergence  constraints.  To 
avoid  imposing  these  constraints  directly,  which  is  difficult,  it  has  been  suggested 
that  they  need  to  be  enforced  only  on  the  boundary  of  the  domain  ([5],  [7] ).  This 
has  some  significant  advantages.  For  instance,  the  vector  Helmholtz  equation 
with  a  divergence  constraint  on  the  field  can  then  be  solved  using  standard  finite 
element  spaces  instead  of  more  complex  spaces  of  edge  elements.  Furthermore, 
as  we  will  show  below,  the  divergence  boundary  condition  can  be  treated  as  a 
natural  boundary  condition. 

It  is  essential  to  know  when  the  interior  and  boundary  formulations  of  diver¬ 
gence  constraints  give  the  same  result  and  that  is  the  main  aim  of  this  work.  We 
shall  show  that  the  formulations  are  equivalent  if  and  only  if  the  scalar  Poisson 
equation  with  any  smooth  right  hand  side  and  Diricblet  boundary  conditions 
has  a  solution  in  When  this  is  not  true,  a  simple  minded  apphcation  of 
divergence  boundary  conditions  as  in  [5]  and  [7]  will  normally  give  incorrect 
solutions  with  a  nonzero  divergence. 

Our  results  for  systems  with  the  interior  divergence  constraint  are  obtained 
by  formulating  the  problem  in  the  space  V  defined  in  the  next  section.  This 
approach  avoids  enforcing  the  divergence  constraint  explicitly  —  it  is  automati¬ 
cally  satisfied.  We  provide  a  proof  of  this  fact  for  interior  problems;  for  exterior 
domains  see  [4],  Unfortunately,  is  not  dense  in  U  in  general.  To  use  stan¬ 
dard  finite  elements  we  impose  the  divergence  condition  on  the  boundary 
and  formulate  the  problem  in  the  space  c  defined  below.  We  discuss 
the  relationship  between  the  V  and  Hq^  formulations  in  some  detail. 

To  avoid  relatively  unimportant  technical  issues  we  will  give  proofe  of  these 
results  for  the  vector  Poisson  equation  subject  to  tangential  boundary  conditions 
and  a  divergence  constraint.  In  the  end  we  will  show  how  they  may  be  extended 
to  other  problems  including  the  vector  Helmholtz  case  of  electromagnetics. 

The  next  section  gives  the  strong  forms  of  the  governing  equations  and  the 
reasoning  behind  the  use  of  divergence  boundary  conditions.  It  also  contains  a 
formula  for  the  boundary  divergence,  which  is  used  to  reformulate  the  divergence 
boundary  condition  as  a  natural  boundary  condition  in  setting.  Section  3 
sets  up  a  weak  form  for  the  equation  with  the  interior  divergence  constraint 
and  proves  the  coercivity  for  the  weak  form,  which  follows  from  a  compact 
embedding  result  for  vector  fields.  Sections  4  and  5  prove  coercivity  for  the  weak 
form  which  uses  divergence  boundary  conditions.  This  time,  coercivity  foOows 
from  a  close  and  somewhat  surprisnig  connection  with  the  previous  weak  form. 
This  result  is  based  on  a  trace  theorem  proved  in  section  5. 

The  coercivity  results  of  sections  3-5  are  used  in  section  6  to  prove  equiv¬ 
alence  of  the  strong  and  weak  formulations  (Theorem  2).  We  also  discuss  in 
what  sense  the  boundary  divergence  condition  is  satisfied  and  prove  a  formula 
for  calculating  the  divergence  on  the  boundary  (Lemma  10  in  section  6) .  Section 
7  contains  a  discussion  of  when  one  can  say  a  priori  that  the  interior  and  bound¬ 
ary  divergence  formulations  are  equivalent,  and  when  they  are  not  equivalent 
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(i.e.  may  have  different  eolutions  for  the  same  data) .  It  also  contains  an  exam¬ 
ple  showing  noneqnivalence  of  the  two  formulations  for  nonconvex  polyhedra. 
In  section  8  we  give  the  modifications  to  include  the  Helmholtz  equation  and 
briefly  mention  some  other  extensions.  Section  9  mentions  some  imphcations 
for  numerical  approximations  of  the  problem. 


2  Formulation  of  the  problems 


Assume  that  H  is  aji  open  bounded  subset  of  ,  JV  =  2  or  3,  with  a  connected, 
Lipchitz  hmmdary  F. 

We  will  consider  the  relationship  between  the  problems 


and 


'-Au  = 

f 

in 

<  div  u  = 

0 

in  f], 

u  X  n  = 

=  0 

on  F 

'-Au  = 

f 

in  n, 

^  div  u  = 

0 

on  F, 

u  X  n  = 

:0 

on  F, 

(1) 


(2) 


We  assume  that  f  6  and  the  compatibility  condition  divf  =  0  is 

satisfied.  Superficially  it  seems  that  the  two  problems  are  equivalent.  Indeed, 
if  u  is  a  solution  of  (1),  it  is  clearly  a  solution  of  (2).  Conversely,  if  u  is  a 
solution  of  (2),  then  talcing  the  divergence  of  the  first  equation  of  (2)  we  get 
that  A{div  u)  =  0,  and  since  div  u  =  0  on  F,  one  might  argue  that  div  u  =  0 
in  n,  so  u  is  a  solution  of  (1).  This  is  basically  the  reasoning  in  [7].  But  it 
might  not  always  be  true.  The  difficulty  is  caused  by  insufficient  regularity 
of  the  solution.  Laplace’s  equation  with  Dirichlet  boundary  conditions  has  a 
unique  solution  in  (Q),  but  in  general  we  cannot  assume  that  div  u  6 
It  may  be  only  in  L^(fl)  and  the  homogeneous  equation  for  div  u  may  have  a 
nonzero  solution  in  L^(H)  (this  happens,  for  example,  in  the  case  of  domains 
with  reentrant  corners). 

In  order  to  discuss  the  the  relationship  between  (1)  and  (2)  we  first  have  to 
make  explicit  in  what  function  spaces  we  look  for  the  solutions.  It  turns  out 
that  the  natural  space  for  (1)  is 

V  =  {ue  L^(n)^  I  divu  6  L^(n),  curlu  e  L^(n)^  ,  u  x  n|r  =  0}, 


where  N'  =  3  H  N  =  S  and  N'  =  I  in  N  =  2.  The  standard  scalar  product  in 
y  is 

(u,  v)y  =  /  u  •  V  +  curlu  •  curl V  +  div udivv. 

Jn 

On  the  other  hand,  the  natural  space  for  (2)  is 


HU^)  =  {ueH\nr  \  MX  n\r=(i} 
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with  the  scalar  product 


(u,v)i  =  /  u  -  ■ 
Jq 


Vu  ■  Vv. 


We  will  show  that  whOe  problem  (1)  has  a  unique  solution  in  it  may  not 
be  solvable  in  i/Q^{n),  and  while  problem  (2)  has  a  unique  solution  in 
the  solution  may  not  he  unique  in  V.  We  will  also  discuss  in  what  sense  the 
divergence  boundary  condition  in  problem  (2)  is  satisfied. 

It  is  shown  in  section  6  (under  the  additional  assumptions  that  T  is  piecewise 
and  the  jump  condition  given  in  (6)  is  satisfied)  that  problem  (2)  (in 
Hq^((\))  is  equivalent  to:  find  u  G 

'-Au  =  f  in  fl, 

^  ^  .  n  _l_  .  n  =  0  on  r,  (3) 

an 

u  X  n  =  0  on  r, 

where  k  is  the  curvature  of  the  boundary  (k  =  div  n),  which  is  defined  almost 
everywhere  on  F.  This  follows  from  the  result  (see  Lemma  10  for  the  precise 
formulation) 

divu|r  =  ^  '  n|r  +  KU  •  n|r. 

It  is  convenient  to  begin  with  the  weak  formulations  of  problems  (1)  and  (3). 
The  next  three  sections  deal  with  the  coercivity  of  the  corresponding  weak  forms. 
Following  this  we  show  the  equivalence  of  the  strong  and  weak  problems  and 
then  examine  the  relationship  of  problems  (1)  and  (3),  and  show  the  equivalence 
of  problems  (2)  and  (3). 


3  Weak  Formulation  in  V 

Let  n  and  f  be  as  above.  To  derive  the  weak  formulation  of  (1)  in  V  we  substitute 
curl  curl  u  for  -Au  in  (1),  multiply  the  first  equation  by  v  (belonging  to  V), 
the  second  by  div  v  and  integrate  by  parts  using  the  boundary  condition.  For 
u,  V  G  y  define 

a  (u,  v)  =  /  curl  u  ■  curl  v  +  div  u  div  v.  (4) 

Jn 

Then  the  statement  of  the  weak  problem  is:  find  u  eV  such  that 

a(u,  v)  =  /  f  •  V  Vv  G  y.  (5) 

in 

The  only  hypothesis  of  the  Lax-Milgram  lemma  which  is  nontrivial  to  check 
is  the  coercivity  of  a  on  V.  It  follows  from  the  next  compact  embedding  theorem. 
The  embedding  is  actually  a  corollary  of  the  regularity  result  in  [1],  where  it  is 
proved  that  V  C  We  have  provided  a  new  and  concise  proof. 
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Theorem  1  V  is  compactly  embedded  in  . 

The  proof  makes  use  of  Muiat^s  div-curl  lemma  (see  [8]). 

Lemma  1  (div-curl  lemma).  Let  U  be  an  open  subset  ofK^.  Let  Vn  and 
Wn  be  weakly  convergent  sequences  in  L^{U)^  with  the  limits  v  and  w  corre¬ 
spondingly.  Assume  that  |divv^}J^i  and  {curlwn}^i  He  in  compact  subsets 
of  [U)  and  respectively.  Then  for  every  (j>  E:'D{U), 

/  (j>Vn  ■  Wn  — /  <j>v  •  w  as  n  oo. 

Ju  Ju 

Proof  of  Theorem  1.  Let  be  a  weakly  convergent  sequence  in  V  with  the 
limit  u.  We  want  to  show  that  converges  strongly  in  L'^{fl).  The  idea  is 
to  extend  outside  fi  in  two  ways:  Vn  and  Wn  will  be  the  extensions  with 
“good”  divergence  and  curl  respectively.  We  construct  the  extensions  so  that 
supp  Vn  n  supp  Wn  c  fl  and  then  use  div-cuxl  lemma  for  these  sequences. 

Let  U  be  an  open  hall  containing  H.  We  start  with  the  construction  of  Vn 
(the  extension  with  divvn  in  a  compact  set  of  H  ^(U)).  For  each  n  6  N  we 
define  gn  6  H^{U\  0)  by 


'Agn  =  0  mU\n, 

<  ^  =  Un  •  n  on  an, 
an 

^5n  =  0  ondU. 


Now  define  by 


and  Wn  by 


Then 


in  n 

mU\n 


Wn 


Un  in  n 
0  in  C/  \  n. 


divVn 


divu^t  in  n, 

0  in  U\n 


and 

^  J  curl  Un  in 

curlwn=<  .  _ 

[o  mU\n 

so  {divvn}^i  {curlwn}JJLi  bounded  subsets  of  L^(C/)  and  , 

respectively  (hence  lie  in  compact  sets  of  and  Moreover, 

Vn  V  in  L^(C/)^  and  Wn  — ^  w  in  Z/^{C/)^,  where  v  and  w  are  similar 
extensions  of  u. 
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Now  choose  (j>  e  VlU)  such  that  <j>  =  l  inQ.  Using  the  div-ciirl  lemma  we 
get  (note  that  =  0  in  U  \  0) 


This  together  with  ^  u  in  shows  that  — >  u  in 

We  will  use  Peetre’s  lemma  (see  [2],  for  example)  in  the  foOowing  form 

Lemma  2  Let  E,  E\  and  E2  he  Banach  spaces.  Let  Ai  and  A2  he  continuous 
linear  operators  from  E  to  Ei  and  from  E  to  E2  respectively.  Assume  that  there 
exists  C  >  0  such  that 

II'i^IIe  <  C  (||^ii/||ei  +  ||^2w||£;2)  Vw.  6  E. 

Assume  also  that  Ker^i  =  {0}  and  that  A2  is  compact.  Then  there  is  Ci  >  0 
such  that 

II^IU  ^  Ml ^^'11  El  Vw  e  E. 

The  coercivity  result  is  then 
Lemma  3  There  is  C  >  0,  independent  of  v,  such  that 
«'(v,  v)  >  ciiviiy  Vv  6  y 

Proof  We  use  Lemma  2  with  E  =  Ei  =  x  L^(n)^  ,  E2  =  1/^(0)^ 
and  ^lU  =  (div  u,  curlu),  A2n  =  u.  By  Lemma  1  the  operator  A2  is  compact. 
We  must  show  that  Ker  Ai  =  {0}.  But 

u  6  Ker  Ai  ^  div  u  =  0  and  curl  u  —  0  in  fl,  and  u  x  n  =  0  on  F. 

This  implies  that  u  =  0  (by  using  the  gradient  potential,  for  example).  So  by 
Lemma  2  the  result  follows. 

The  Lax-MUgram  lemma  now  implies  that  (5)  has  a  unique,  stable  solution 
in  V. 


4  Weak  Formulation  in 

For  deriving  the  weak  form  of  (3)  we  shall  make  additional  assumptions  about 
the  smoothness  of  the  boundary.  In  the  following  11  is  a  hoimded  subset  of  , 
N  =  2  or  3,  with  a  connected,  Lipchitz  and  piecewise  boundary  F,  i.e. 

F  =  uy^iFj  with  Fj  n  fi=  0  for  2  7^  j,  and  F^,  j  =  1, . . .  ,n  are  of  class 
Let  7  denote  the  set  of  ‘^edges  and  comers’’,  i.e. 

n 

7=  U  (rinr,). 

■*  J=1 

i^j 
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We  assume  that  the  jump  of  the  normal  on  7  is  bounded  below,  i.e.  there  is 
5  >  0  such  that  for  all  x  €  7  with  x  6  n  i  ^  j  we  have 

n^(x)  •  n^  (x)  <1-5,  (6) 


where  n*  (x)  and  n’  (x)  are  the  limits  of  the  unit  outer  normals  when  approaching 
X  from  Ti  and  F^- ,  correspondingly. 

This  condition  is  satisfied  for  all  polygons  and  polyhedra.  Also  included  are 
polygons  and  polyhedra  with  curved  sides.  It  excludes  three-dimensional  bodies 
which  have  points  like  the  tip  of  a  cone,  and  bodies  with  edges  which  “flatten 
out”,  i.e.  the  angle  between  the  faces  gets  arbitrarOy  close  to  tt.  This  condition 
is  needed  for  proving  some  results  about  trace  operators  in  the  next  section, 
and  imphcitly  in  the  proof  of  Lemma  4,  which  gives  the  relation  between  the 
two  weak  forms. 

To  obtain  a  suitable  weak  form  for  (3)  we  proceed  formally,  multiplying  the 
equation  by  a  test  function  v  6  and  integrating  by  parts  as  usual.  We 

get 

f  Vu-Vv-  f  [  f  •  V. 

Jn  Jr  Jn 

Rewriting  the  boundary  term  as  a  sum  over  the  smooth  boundary  pieces  and 
using  the  boundary  conditions  on  u  and  v  gives  the  weak  problem:  find  u  € 
such  that 


f  Vu-Vv-h^/  «:u-v=  ^f-v  VvGiJot(^)- 
Ja  Jrj  Jn 


(7) 


Since  «  6  L'^(Fj),  by  the  usual  trace  theorem  the  boundary  term  is  well  de¬ 
fined.  Note  that,  in  this  framework,  the  divergence  boundary  condition  (or  the 
equivalent  one  in  problem  (3))  is  natural. 

Proving  coercivity  of  the  weak  form  (7)  is  a  nontrivial  matter  because  the 
curvature  can  be  of  either  sign.  The  result  will  follow  from  the  next  lemma 
whose  proof  is  deferred  to  the  foUowing  section. 


Lemma  4  Assume  that  H  satisfies  ike  assumptions  made  in  the  beginning  of 
the  secMon.  For  u,v  6  we  have 

a(u,  v)  =  /  Vu  » Vv  +  f  KU  ■  V,  (8) 

Ja  ^^Jtj 


where  F^-,  j  =  1, . . .  ,  n  are  the  smooth  pieces  ofT,  and  a  is  the  same  as 

in  (4). 

This  means  that  the  weak  problem  (7)  is  in  fact  equivalent  to:  Find  u  e 
Buch  that 

a(u,v)=  /f-v  (9) 

Jn 
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in  which  the  equation  ie  the  same  as  in  (5)  but  the  spaces  axe  different. 

In  fact,  one  can  get  (9)  from  problem  (1),  but  they  are  not  equivalent: 
in  general  we  cannot  assiime  that  the  solution  of  (1)  belongs  to  and 

we  cannot  prove  later  that  the  solution  of  (9)  satisfies  divu  =  0.  We  can 
also  formally  get  (9)  from  (2)  by  substituting  —  Au  by  curl  curl  u  —  Vdivu, 
multiplying  by  v  and  integrating  by  parts  using  the  boundary  conditions,  but 
this  is  not  easy  to  justify  rigorously  (for  the  solution  of  (2),  curl  curl  u  and 
Vdiv  u  may  not  be  in  L^). 

The  coercivity  result  for  (9)  is: 

Lemma  5  There  is  (7  >  0  such  that 

a(u,v)>C||u||i  Vu€F^,(n). 

Proof.  We  use  Lemma  2  with  E  =  Ei  =  L^(Q)  x  ,  E2  = 

L^(n)^  X  L^(r)^,  Aiu  =  (divu,curlu)  and  A2n  =  (u,u|r).  The  operator 
A2  ■  E  ^  E2  ^  compact  because  it  is  bounded  as  an  operator  from  E  to 
X  and  the  latter  is  compactly  embedded  into  E2.  We  already 

showed  (in  the  proof  of  Lemma  3)  that  Ker  Ai  =  {0}.  The  result  follows. 

Applying  the  Lax-MUgram  lemma  to  the  weak  form  (9)  (or,  equivalently,  to 
(9))  and  using  the  coercivity  result  above  shows  the  existence  of  a  unique  stable 
solution  to  (9). 

Note  that  to  prove  the  coercivity  of  a.  in  we  do  not  need  the  additional 

smoothness  assumptions  made  in  the  beginning  of  this  section:  they  are  needed 
to  mahe  sense  of  the  term  containing  the  curvature  in  (7)  and  to  show  later  that 
the  divergence  boundary  condition  is  satisfied  for  the  solution  of  problem  (9). 
They  are  not  needed  for  uniqueness  of  solution  of  (9). 


5  Proofs  for  previous  section 


This  section  will  present  the  proof  of  Lemma  4.  We  will  assume  throughout  the 
section  that  the  assumptions  made  in  the  beginning  of  the  previous  section  are 
satisfied. 

First  we  need  some  preliminary  results.  We  begin  by  recalling  some  facts 
about  the  trace  spaces.  First,  for  Lipchitz  F  c  the  following  is  an  equivalent 
norm  on  (see  e.g.  [3]): 

II,AI|2  -  ll^l|2  A-  f  f  rfC 


We  will  also  need  the  spaces  Fo^^(Fj),  j  =  1, . . .  ,n,  which  consist  of  all  func¬ 
tions  in  iJ^/^(Fj)  whose  extension  by  zero  to  F  belongs  to  i7^/^(r).  The  norm 
of  a  function  in  is  the  norm  of  its  extension  by  zero  in  i7^/^(r),  i.e. 


m'- 


+2/ 

Mr,  M  1^  “ 


jY  dSjf^dSy. 
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We  wUl  use  the  same  notation  for  a  function  in  HQQ^{rj)  and  its  extension  by 
zero  to  r. 

We  will  also  need  the  following  result: 


Lemma  6  The  trace  operator  v  v-n|rj  is  a  continuous  linear  operator  from 
onto 

Proof.  We  use  the  fact  that  the  usual  trace  operator  v  v|r  is  a  continuous 
linear  operator  from  (H)^  onto  (F)  ^ .  First  we  will  show  that  if  v  x  n  =  0 
and  the  jump  condition  (6)  is  satisfied,  then  ll’^•^lli^V2(r^.)  <  ^l|v||//V2(r)-  Let 
^  =  V  •  n  on  r^*  and  ^  =  0  on  F  \  F^.  Clearly  ||^|i/fV2(r^.)  <  ||v|l/iri/2(r^.),  so  we 
have  to  show  only  that 

/\r,  /r, 

Because  of  the  jump  condition  we  can  find  e  >  0  such  that 

X  €  r^,  y  6  r  \  Fj,  |x  -  y|  <  e  =>■  n(x)  •  n(y)  <  1  “  g  • 


For  X  €  F^-,  y  6  F  \  F^  such  that  |x  —  y|  <  e  we  have 
|v(x)  -  v(y)|^  =  |v(x)p  +  |v(y)|^  -  2(v(x)  ■  n(x))(v(y)  ■  n(y))(n(x)  •  ii(y)) 

>  |«/i(x)p  +  |v(y)p  -  2|^(x)|  |v(y)|  (l  -  0 

>  ^1  -  j  1^5  -  . 


6^ 

Since  S  — ^  ^  (note  that  0  <  5  <  2,  otherwise  the  jump  condition  could  not 
be  satisfied),  we  have 
I^(x)i2 

/r\rj  Jtj 


f  f  l^(x 
Jr\ri  Jvj  “ : 

=/  / 


jY  dSxdSy 


AN 


dS^dS. 


lSy+  f  ( 

JT\T.j  Jr. 


l^(x)P 


•jnB4y)  |x  -  y| 

v{x)-v{y)|2 


/r\r,-  Jrj\B4y)  |x-y| 

I  [ 

^  yr\rj-  Jrj  ^  “  T  J^\^3 


_  ^\N 


■yl 


N 


Consequently,  the  operator  v  v  •  nlr^.  is  a  continuous  linear  operator  from 

HM  to 

Let  <j>  €  HQQ^{rj)  be  given.  Ebctend  it  by  0  to  F.  Then  and 

we  can  find  v  G  77^(0)^  such  that  v|r  =  and  |lv||i  <  C'||^/>n||j:^i/2(r)  ^ 


8 


Cill^ll  „i/2.p  ..  Clearly  v  6  so  the  trace  operator  is  onto. 

^00 

The  next  lemma  gives  a  formula  for  calculating  the  divergence  on  the  bound¬ 
ary 


Lemma  7  For  smooth  u  the  restnction  of  div  u  to  Vj  satisfies 


divu]r^.  =  (  divrUr  +  KU  ■  n  +  —  •  n  I 


where  is  the  tangential  component  of  u  and  divp  is  the  divergence  of  a  tan¬ 
gential  vector  field  in  the  tangential  coordinate  system. 

Proof  To  calculate  div  u  at  some  xq  G  let  us  fix  a  coordinate  system  i,  j, 
k  (if  iV  =  2  then  omit  j),  where  i,  j  are  tangent  to  F  at  xq  and  k  is  the  normal  at 
xo  (we  can  do  this  for  every  xq  in  the  interior  of  Tj).  Denote  =  u  ~  (u  ■  n)n. 
Then 

divu(xo)  =  ^(xo)  •  1  +  ^(xo)  •  J  +  ^(^o)  •  k 

dUr ,  .  .  dur.  .  .  9[(u-n)n]  .  5[(u-n)n]  . 


+  |^(5Co)  ■  n(xo) 


=  (di 


V 


divrUr  +  Ku  •  n  + 


•”) 


where 


divrUT(xo)  =  ^(xn)  •  i  +  ^(^o)  •  j. 


Now  we  can  prove  lemma  4: 

Proof  of  Lemma  4.  For  v  G  aiid  smooth  u  (we  do  not  require 

u  X  n  =  0  yet,  because  it  does  not  seem  straightforward  to  show  the  density  of 
smooth  functions  with  zero  tangential  trace  in  HqX^)) 

/  curl u  ■  curl V  +  div u div  V  =  /  curlcurlu- v  —  Vdivu  •  v  +  /(divu)v  n 

Jn  Jr 

=  —  Au  •  V  4-  y  (div  u)v  ■  n 


—  [  Vu  -Vv+  [  I  divu  —  ^  •  n  I  V  •  n. 

Ju  Jr\  dn  J 


By  Lemma  7 


/  curl u  ■  curl V  +  div u divv  =  /  Vu-Vv  +  Y^  /  (divrUT  +  ku  •  n)  v  •  n. 

Jn  Jn  (10) 
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For  any  <j>  e  (extended  by  zero  to  F)  we  can  find  v  €  such 

that  vjr  =  and  ||v|li  <  C\\<j>\\^ij2^^^y  EVom  (10)  we  get 


I 


(divrUr  +  «:u  •  n)  ^ 


6  Vu  smooth. 


By  density  of  smooth  functions  in  H^{D)  and  Hahn-Banach  theorem  we  can 
continue  the  mapping  u  (diYru^  +  «;u  ■  n)|rj  uniquely  to  a  continuous  linear 

mapping  from  to  (i?oo^(Fj))'.  Since  for  any  v  G  v  •  n  G 

(by  Lemma  6),  (10)  holds  for  all  u  G  and  v  G  (f^^ 

integrals  over  F^  have  a  meaning  as  a  duality  pairing  between  (i?oo^(Fj))'  and 

To  finish  the  proof  we  need  to  show  only  that  if  u  G  H^f{u)y  then  divrti^  =  0 
on  F^  ,  For  this  first  note  that  since  u  «u  •  n|r^.  is  hnear  and  continuous  from 
(H)^  to  c  (Fo^^(F^))',  the  operator  u  divpu,-  Ir^.  is  also  continuous 

from  to  and  then  use  the  fact  that 

(divru,,  <l>)r,  =  -  /  u,  •  €  Cl  {Tj) 

Jrj 

and  the  density  of  Cl  (Fj)  in 


6  Equivalence  of  the  Strong  and  Weak  Forms 

In  this  section  we  will  prove  that  (1)  in  V  is  equivalent  to  the  weak  problem 
(5),  and  that  (3)  in  is  equivalent  to  the  weak  problem  (7).  It  follows 

that  problem  (1)  is  uniquely  solvable  in  V  and  problem  (3)  is  uniquely  solvable 
in  i7o£(n).  We  vdll  also  show  that  problems  (2)  and  (3)  are  equivalent  with  the 
divergence  boundary  condition  having  a  meaning  in  the  sense  of  traces. 

Let  us  first  deal  with  problem  (1). 

Lemma  8  Let  fl  and  f  he  as  in  Section  2.  Then  the  problem  of  finding  ueV 

satisfying  (1)  is  equivalent  to  solving  the  weak  problem  (5). 

Proof  Let  u  G  y  he  a  solution  of  (1).  Then  -Au  =  curl  curl  u  G  L‘^{Q)) 
and  so  the  formal  calculations  leading  to  (5)  can  be  rigorously  justified. 

Conversely,  if  u  G  V  satisfies  (5),  then  by  using  v  G  P{n)^  and  transferring 
all  derivatives  to  v  we  get  —  Au  =  f  in  the  sense  of  distributions.  For  any 
(j>  G  'D{Q)  we  can  find  g  e  such  that  Ag  =  Now  use  v  =  in  the 

weak  form  (note  that  Vg  eV).  We  get 

f  (div  u)^  =  [  f  ■  =  -  /  (div  f)^  +  f  (f  ■  n)p  =  0, 

Jn  Jfi  Jn  Jr 
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80  div  u  =  0  in  the  sense  of  distributions.  Consequently  the  problem  (1)  (in  V) 
and  the  weak  form  (5)  are  equivalent. 

Corollary  1  The  problem  (1)  has  a  unique  solution  Uv  in  V. 

In  the  proof  of  the  theorem  we  used  the  fact  that  Vg  eV.  If  we  were  dealing 
with  the  weak  form  in  then  Vg  may  not  have  been  a  legitimate  test 

function  (in  the  case  g  ^  which  would  make  it  impossible  to  prove  that 

div  u  =  0.  This  is  the  difference  between  the  weak  problems  (5)  and  (9),  even 
though  they  appear  very  simOar. 

The  corresponding  proof  for  the  problem  (3)  is  a  little  more  difSlcult,  bnt 
foOows  the  same  pattern. 

Lemma  9  Let  11  and  f  he  as  in  Section  2.  Assume  in  addition  that  the  boundary 
r  is  piecewise  and  the  jump  condition  (6)  is  satisfied.  Then  the  problem  of 

finding  u  6  satisfying  (S)  is  equivalent  to  solving  the  weak  problem  (7). 

Proof.  We  first  have  to  show  that  the  boundary  condition 

—  •n4-/tu-n=0  onF 
on 

has  a  meaning  for  u  G  with  Au  6  Since  Vu  G  and 

divVu  G  the  gradient  has  a  normal  trace  on  the  boundary  —  G 

H  ^/^(F)  (see  [2]  for  example).  For  <j>  G  we  have  G 

(using  the  extension  of  <j>  by  zero)  and  therefore  the  boundary  condition  may  be 
interpreted  as 

+  «u,<^n)r  =  0  e  n  (11) 

.7  =  1 

or,  equivalently, 

+  «U,  v)r  =  0  Vv  €  (12) 

where  {■,  •)r  is  the  duality  pairing  between  7f"^/^(F)  and 

Let  u  G  ^  solution  of  (3).  Using  (12)  we  can  justify  to  formal  cal¬ 

culations  leading  to  the  weak  form  (7)  (the  boundary  integrals  in  the  weak  form 

should  be  understood  as  duality  pairings  between  (i^oo^(^7))^  ^oo^(^i))- 

Conversely,  if  u  G  (7),  then  by  using  v  G  P(n)^  we  get  as 

before  -  Au  =  f  in  the  sense  of  distributions.  For  any  v  G  LTd*  (0),  integrating 
by  parts  in  the  weak  form  (note  that  Au  G  1/^(0)^)  we  get 

hence  (12)  is  satisfied.  Consequently  the  problem  (3)  (in  i7d^  (fi))  and  the  weak 
form  (7)  are  equivalent. 
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Corollary  2  The  problem  (3)  has  a  unique  solution  u  in 

Now  we  will  show  that  problems  (2)  and  (3)  are  equivalent.  The  only  dif¬ 
ference  in  these  two  problems  is  in  the  boundary  condition,  so  we  have  to  show 
that  the  divergence  boundary  condition  makes  sense  and  is  equivalent  to  the 
boundary  conditon  in  (3).  This  is  done  in  the  following  lemma. 

Lemma  10  Assume  H  satisfies  ike  conditions  of  Lemma  9.  For  u  6 
with  Au  6  mar,  and  any  i  6  {1, . . .  the  trace  o/div  u  on  Tj  exists  and 
belongs  to  {HQQ^{rj)y ,  and 

.  .  du  .  . 

^  ■  ^lr,-  ♦ 

Proof.  This  follows  from  Lemma  (7)  which  gives  the  formula  for  smooth  u 
without  the  condition  ux  n  =  0.  By  density  this  can  be  extended  ton  e 
with  Au  €  .  Then  use  the  fact  that  for  u  G  Hq^{Q)  we  have  div  =  0 

(see  the  end  of  proof  of  Lemma  4). 

Corollary  3  Problems  (2)  and  (3)  are  equivalent. 

Let  us  summarize  the  results  in  a  theorem: 

Theorem  2  Assume  that  fl  is  a  bounded  subset  of  vdth  a  connected,  Lip- 
schitz  boundary.  Let  f  G  mar  with  div  f  =  0  given.  Then 

a)  The  problem  (1)  inV  is  equivalent  to  the  weak  form  (5)  and  has  a  unique 
solution  ny  G  V. 

b)  If  in  addition  the  boundary  is  pieceuhse  and  satisfies  the  jump  cori- 
dition  (6),  then  the  problems  (2)  and  (3)  in  H^^{Q)  are  both  equivalent  to  weak 
problems  (7)  and  (9),  and  have  a  unique  solution  n^  G 


7  Relationship  of  the  V  and  -f^ot(O)  Formulations 

Let  uy  G  y  be  the  solution  of  (5)  (or  the  strong  form  (1  with  the  interior 
divergence  condition)  and  let  6  be  the  solution  of  (9)  (which  under 

the  additional  smoothness  conditions  is  equivalent  to  problems  (2)  and  (3)  with 
the  divergence  boundary  condition).  In  this  section  we  will  discuss  when  the 
two  solutions  are  the  same.  The  following  is  clearly  true. 

Lemma  11  Let  n  and  f  be  as  in  Section  2.  Assume  in  addition  that  the  jump 
condition  (6)  is  satisfied.  Then  the  following  are  equivalent 
a)  Uy  =  Uh; 
h)  div  Uff  =  0  in  fl; 

c)  Uv  e 

Let  us  now  examine  when  uy  =  Uh.  First,  note  that  lemma  5  implies 
that  HQ^{Q)  is  closed  in  V.  We  have  two  possibilities:  either  V  = 
or  y  ^  i?^^(n).  In  the  first  case  the  solutions  are  obviously  the  same  (the 
corresponding  weak  forms  are  exactly  the  same).  The  following  lemma  gives 
the  necessary  and  sufficient  conditions  for  the  equality  of  the  spaces. 
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Lemma  12  Assume  that  fl  is  a  bounded  subset  with  a  connected  Lip- 

schitz  boundary.  The  spaces  V  and  are  equal  if  and  only  if  the  homo¬ 

geneous  DiHchlet  problem  for  the  scalar  Poisson  equation  with  any  right  hand 
side  in  V{Q)  has  a  solution  in 

Proof  Suppose  that  the  regularity  condition  holds.  For  any  given  w  G  V" 
we  can  solve  the  weak  problem:  find  u  e  (H)  such  that 

a(u,  v)  =  a(w,  v)  Vv  €  Hot (^)  (1®) 

(Le.  u  is  a  projection  of  w  onto  corresponding  to  the  scalar  product 

a(-,  •)  ).  As  in  the  proof  of  Lemma  8  we  can  ^ow  that  A(u  -  w)  =  0,  and 
using  the  i7^-regularity  of  the  solution  of  the  Dixichlet  problem  for  the  Poisson 
equation  we  also  get  div  (u  —  w)  =  0.  Since  (u— w)  x  n|r  =  0,  by  the  uniqueness 
of  solution  of  (1)  we  have  u  =  w,  Le.  w  6  Consequently  V  =  HQ\{n). 

Now  suppose  that  there  is  p  6  such  that  Ag  €  P(n),  but  g  ^ 

Then  Vg  €  F,  but  Vg  ^  soV  ^ 

Note  that  the  lemma  remains  true  if  one  substitutes  for  'D{Q). 

Corollary  4  J/L!  is  convex  or  has  a  boundary  then  V  — 

This  follows  from  the  classical  regularity  theory  about  the  smoothness  of  the 
solution  of  the  Dirichlet  problem  for  the  Poisson  equation.  It  is  actually  enough 
to  assume  that  H  is  locally  convex  near  the  points  where  the  boundary  is  not 
(i.e.  these  points  have  a  neighborhood  whose  intersection  with  His  convex). 

The  question  of  the  equahty  of  these  spaces  (and  analogous  ones  with  a 
boundary  condition  for  the  norrnal  trace)  has  been  partly  answered  by  different 
authors,  but  we  are  unaware  of  results  as  sharp  as  these  presented  here.  In  the 
classical  book  [2]  it  is  proved  that  V  =  if  Q  is  a  convex  polygon  or  has 

a  boundary.  Necessary  and  sufficient  conditions  for  V  —  to  hold  in 

are  given  in  [6],  namely  that  the  scalar  Poisson  equation  with  L^(n)  right 
hand  side  must  have  a  solution  in  77^(0)  both  with  Dirichlet  and  Neumann 
condition.  As  we  saw,  the  regularity  of  solution  of  the  Neumann  problem  is 
not  needed,  and  in  it  actually  follows  from  the  regularity  of  solution  of  the 
Dirichlet  problem. 

In  the  other  case,  when  ^  identical,  since  Hq^{Q)  is 

closed  in  V,  there  is  f  6  L^(n)  with  div  f  =  0  for  which  the  solution  of  (1)  is  in 
y,  but  not  in  means  that  the  interior  and  boundary  divergence 

formulations  axe  not  equivalent  and  give  different  solutions  for  this  f .  We  wiU 
show  how  to  construct  f  for  which  uh  whenever  V  7^ 

Lemma  13  IfV^  than  there  is  f  6  L^{Q)  with  div  f  =  0  such  that  the 

corresponding  solution  of  (1)  is  not  in 

Proof.  If  y  7^  (0)  then  we  can  choose  a  nonzero  w  6  y  such  that 

a(w,  v)  =  0  Vv  €  (14) 
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First  note  that  curl  w  7^  0:  indeed,  if  we  had  curl  w  =  0  then 
div  w  •  div  V  =  0  Vv  €  ^ 


i 


and  hence  (by  using  v  6  2^(0) Vdivw  =  0,  01  divw  =coiiBt.  It  follows 
from  (15)  that  div  v  =  0  for  all  v  6  {^)i  which  is  a  contradiction. 

Now  let  f  G  y  be  a  solution  of 


a(f,v)=  f 


curl  w  ■  curl  v  Vv  G  V. 


(16) 


Note  that  f  7^  0  and  divf  =  0.  We  claim  that  for  this  f  problem  (1)  has  a 
solution  uy  G  y  which  is  not  in  Indeed,  using  the  weak  problem  (5) 

and  (16)  we  get 

a(u,w)  =  [  curl  u  •  curl  w  —  a(u,  f)  =  f  |f|^  7^  0, 

Jfi  Jn 

which  by  (14)  implies  that  uy  ^ 

We  can  also  give  a  more  tangible  example  of  a  solution  of  (1)  which  does 
not  belong  to  i7^^(n).  Suppose  n  C  is  such  that  in  a  neighborhood  of  some 
xo  G  r  its  boundary  consists  of  2  planes  with  the  interior  angle  a  >  tt.  Choose 
a  cylindrical  coordinate  system  {r,9,z)  with  the  origin  at  xq,  z-axis  along  the 
edge  and  ^  =  0,  0  =  a  on  the  boundary.  Let  ^  G  C'^(n),  depending  on  r  and 
js;  only,  be  such  that  if>=  1  near  the  origin  and  ^  =  0  outside  a  neighborhood  of 
the  origin  where  the  pieces  of  F  are  planar.  Put 

u(x,  y,  z)  =  curl  ^^(r,  z)r'^^°^  cos  ^  , 

where  is  the  unit  vector  in  direction  of  2;-axis.  Then  div  u  =  0  and 


u  X  n 


dn 


z)r^^^  cos  — ^  Bz  =0 


on  r. 


It  foIlowB  that  u  ifl  a  solution  of  (1)  with  f  =  —  Au  G  1/^(0)  and  divf  =  0.  In 
addition,  f  =  0  near  the  origin.  Since  u  ^  and  the  problem  (2)  with  the 

same  f  has  a  solution  u  G  we  must  have  u  7^  u  and  div  u  7^  0. 

This  example  is  essentially  two-dimensional.  To  make  it  work  for  fl  C  R^ 
in  a  neighborhood  of  a  corner  with  the  interior  angle  or  >  tt  we  just  use  <}> 
independent  of  2;.  Note  that  the  third  component  of  u  is  zero,  and  we  get 

f  G  C^(n)2. 

We  will  state  these  results  as  a  theorem: 

Theorem  3  Assume  that  fl  is  a  hounded  subset  of  R^  with  a  Lipschitz  bound¬ 
ary.  The  following  are  equivalent: 

a)  uv  =  u//  for  any  f  G  L^(Q)^  with  divf  =  0; 

c)  the  DiHchlet  problem  for  the  scalar  Poisson  equation  with  any  smooth 
Hght  hand  side  has  a  solution  in 
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8  Generalization  to  Helmholtz  equation 

Most  of  the  argument B  above  remain,  true  in  case  of  Helmholtz  equation.  The 
only  difference  is  in  the  coercivity,  which  now  will  not  be  sufficient  for  using 
Lax-Milgram.  Let  us  sketch  the  results  briefly.  The  analog  of  problem  (1)  is 
now:  find  u  6  F  such  that 

— Au  —  k'^u  =  f  in  f], 

^  divu  =  0  in  n,  (17) 

u  X  n  =  0  on  r. 

We  assume  that  f  6  L^(n)^  and  the  compatibUity  condition  divf  =  0  is  satis¬ 
fied.  To  formulate  the  weak  problem  let  us  define 

6(u,  v)  “  /  curl  u  •  curl  v  +  div u div  v  —  k^n  ■  v  Vu,  v  6  V. 

Jn 

Then  the  weak  problem  corresponding  to  (17)  is  to  find  ueV  such  that 

6(u,v)rr  [  f.v  VvGT^.  (18) 

The  coercivity  result  is  that  &(u,  u)  >  Cl|u||y  —  A:^||u||^2  (if  is  complex, 
we  should  use  the  real  part).  We  cannot  use  Lax-MUgram  directly.  Instead 
we  can  use  the  Fredholm  alternative  to  show  the  existence  and  uniqueness  of 
the  solution.  Define  operator  Ay  :  ^^(0)  L^{Q)  by  the  requirement  that 

Ayg  G  V  is  the  solution  of  a(u,  v)  =  /n  g  •  v  Vv  6  F  (note  that  we  do  not  need 
the  compatibility  condition  div  g  =  0  to  guarantee  the  solvabihty  of  the  weak 
problem).  By  the  Lax-MUgram  Lemma  and  the  compact  embedding  theorem 
(Theorem  1)  Ay  is  a  compact  operator.  Then  the  solution  of  (18)  is  the  solution 
of  u  =  Ay  {k'^\i  +  f),  and  if  1/A:^  is  not  an  eigenvalue  of  Ay  then  (18)  has  a 
unique  solution  ueV.  The  operator  Ay  has  two  kinds  of  eigenfunctions:  ones 
that  are  divergence-free,  and  others  which  are  gradients  of  the  eigenfunctions  of 
the  scalar  laplacian  with  Dirichlet  boundary  conditions.  Assuming  that  1/A:^  is 
not  an  eigenvalue  of  Ay ,  one  can  prove  the  equivalence  of  the  strong  and  weak 
problems  simUarly  to  the  proof  of  Lemma  8  (in  proving  that  the  solution  of  the 
weak  form  satisfies  the  divergence  condition  we  need  solvability  of  the  scalar 
Helmholtz  equation  with  Dirichlet  boundary  conditions,  which  follows  from  the 
fact  that  if  1/A:^  is  not  an  eigenvalue  of  Ay ,  then  it  is  not  an  eigenvalue  of  the 
scalar  laplacian  with  the  Dirichlet  boundary  conditions). 

The  analog  of  problem  (2)  is  to  find  u  g  Ho£(^)  such  that 

— Au  — /c^u  =  f  in  n, 

<  div  u  =  0  on  r,  (19) 

u  X  n  =  0  on  r 

and  the  corresponding  weak  problem  is  to  find  u  G  Hot(^)  such  that 

6(u,  v)  =  /  f  •  V  Vv  G  Ho\(n).  (20) 

Jn 
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One  can  deal  with  the  coercivity  as  above,  using  the  operator  :  L2(0) 
L^{(1)  defined  by:  Ang  €  Fo((n)  is  the  solution  of  «(u,v)  =  /g  vVv  €  /fot(n). 
Equi'^ence  of  the  strong  and  weak  forms  foUows  similarly  to  the  zero  frequency 
case. 

The  question  when  the  two  weak  problems  (18)  and  (20)  give  the  same 
solution  can  he  answered  as  in  section  7.  Most  of  the  results  presented  there  do 
not  depend  on  the  frequency.  The  counterexamples  are  also  easy  to  modify:  if 
uy  ^  18  a  solution  of  (1),  then  it  is  also  a  solution  of  (17)  with  the  right 

hand  side  f  -  k’^uy .  Theorem  3  remains  true  if  we  assume  in  addition  that  1/fc^ 
is  not  an  eigenvalue  of  either  Ay  or  Ah  and  let  uy  and  uh  be  the  solutions  of 
the  weak  problems  (18)  and  (20)  correspondingly. 

Similar  results  can  be  proved  for  time-dependent  problems,  such  as  vec¬ 
tor  wave  equation  and  Maxwell's  equations.  Even  with  divergence-free  initial 
conditions,  using  a  we^  form  in  Hq^{Q)  (similar  to  (20))  may  give  a  solution 
not  satisfying  the  divergence  constraint.  Again,  this  can  happen  only  if  the 
iJ^-regularity  of  the  Dirichlet  problem  for  the  scalar  Poisson  equation  fails. 


9  Conclusions 

We  saw  that  in  some  cases  (e.g.  for  noDconvex  polygons)  the  boundary  diver- 
gence  formulation  is  not  equivalent  to  the  interior  divergence  formulation,  and 
the  weah  form  (9)  or  (20)  (the  weak  forms  in  uniquely  solvable, 

may  not  have  a  divergence-free  solution.  In  this  case  all  approximation  meth¬ 
ods  based  on  these  weak  formulations  (e.g,  the  usual  finite  element  method, 
the  least  squares  method)  will  converge  to  a  spurious  solution  (not  satisfying 
the  divergence  constraint).  Tb  avoid  the  spurious  solution,  one  should  use  the 
weak  form  in  V,  but  this  is  not  straightforward:  one  has  to  use  basis  functions 
which  are  in  F,  but  not  in  i7^t(n)  (e.g.  the  singular  solutions  around  the  re¬ 
entrant  corners),  in  addition  to  the  usual  ones.  The  same  is  true  about  the 
penalty  method,  where  the  weak  form  contains  a  penalty  parameter  in  front 
of  the  terra  with  divergence.  H  the  solution  of  the  original  problem  is  not  in 
then  any  choice  of  the  penalty  parameter  will  result  in  a  spurious  solu¬ 
tion  —  in  fact,  the  smaller  one  makes  the  divergence  (by  choosing  larger  values 
of  the  parametei),  the  bigger  the  error  is  in  the  curl  of  the  solution. 

On  the  other  hand,  when  the  formulations  are  equivalent  (this  is  the  ease 
when  the  solution  of  the  Dirichlet  problem  for  the  scalar  Poisson  equation  is  in 
then  one  can  use  the  weak  forms  in  enables  one  to  use 

simpler  algorithms,  e.g.  finite  element  methods  with  piecewise  linear  test  and 
trial  functions.  The  equivalent  formulation  in  terms  of  gradients  (as  in  (7))  may 
be  especially  useful,  since  it  decouples  the  field  components  inside  the  domain. 
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